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Abstract 



We present a detailed analysis of e' je within the Standard Model, taking into ac- 
5^ , count the strong enhancement through final-state interactions identified in refs. Q 

and The relevant hadronic matrix elements are fixed at leading order in the 
l/Nc expansion, through a matching procedure between the effective short-distance 
Lagrangian and its corresponding low-energy description in Chiral Perturbation The- 
ory. All large logarithms are summed up, both at short and long distances. Two 
different numerical analyses are performed, using either the experimental or the the- 
oretical value for e, with compatible results. We obtain Re (e'/^) = (1.7ib0.9) • 10~^ . 
The error is dominated by the uncertainty in the value of the strange quark mass 
and the estimated corrections from unknown l/A'c-suppressed local contributions. 
A better estimate of the strange quark mass would reduce the uncertainty to about 
30%. The Standard Model prediction agrees with the present experimental world 
average Re (e'/e) = (1.93 ± 0.24) • 10"^ . 



1 Introduction 



The CP-violating ratio e' je constitutes a fundamental test for our understanding of fla- 
vour-changing phenomena within the Standard Model framework. It represents a great 
source of inspiration for physics research and has motivated in recent years a very inter- 
esting scientific controversy, both on the experimental and theoretical sides. 

The experimental status has been clarified recently. The CERN NA48 collaboration 
has announced a preliminary value 

Re(e7e) = (1.40±0.43) ■ 10"^ . (1) 
A larger result was obtained by the Fermilab KTeV collaboration [Q, 

Re (e'/e) = (2.80 ±0.41) ■ 10^^ . (2) 
The present world average 0, 

Re (e'/e) = (1.93 ±0.24) ■ 10^3 , (3) 

provides clear evidence for a non-zero value and, therefore, direct CP violation phenomena. 

The theoretical status is more involved and not very satisfactory. There is no universal 
agreement on the e' je value predicted by the Standard Model, since different groups, using 
different models or approximations, obtain different results [p|-p!^. Nevertheless, it has 
been often claimed that the Standard Model predicts a too small value of e'/e, failing 
to reproduce its experimental world average by at least a factor of two. This claim has 
generated a very intense theoretical activity, searching for new sources of CP violation 
beyond the Standard Model framework [113 . 



It has been pointed out that the theoretical short-distance evaluations of e' je had 
overlooked the important role of final-state interactions (FSI) in ^ tttt decays. Although 
it has been known for more than a decade that the rescattering of the two final pions induces 
a large correction to the isospin-zero decay amplitude, this effect was not taken properly 
into account in the theoretical predictions. From the measured vr-vr phase shifts one can 
easily infer that FSI generate a strong enhancement of the predicted e' j e value, by roughly 
the needed factor of two [1|,0]. A detailed analysis of the corrections induced by FSI has 
been already given in ref. 0, where the low-energy (infrared) physics involved has been 
investigated and the size of the FSI enhancement and the associated uncertainties have 
been quantified. 

In this paper, we present a complete reevaluation of e' je within the Standard Model. 
We will show that with our present understanding of the different inputs, it is possible to 
pin down the prediction of this important parameter with a theoretical accuracy of about 
50%. In order to achieve this goal, one needs to identify the most important corrections 
and find appropriate expansion parameters to perform a perturbative approach with well- 
defined power counting. 

The large- A^^c expansion |1^,|T^, with A'',^ the number of QCD colours, turns out to be a 



very useful tool to organize the calculation. It is a unique non-perturbative approach, with 
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a clear meaning within the usual perturbative expansion in powers of the QCD coupling. 
At leading (non-trivial) order in 1/Nc it is possible to compute all needed ingredients and, 
what is even more important, the matching between short- and long-distance physics can 
be done exactly. Moreover, FSI are zero at leading order in 1/Nc] this allows a clear 
separation of these corrections, avoiding any possible ambiguity or double-counting. 

Since Nq = 3 in the real world, the natural size to be expected for the l/Ai'c-suppressed 
contributions is 30%. Actually, there is a quite compelling phenomenological evidence that 
those corrections are usually smaller. For this to be true, however, one needs to make sure 



that the 1/Nc expansion does not involve large logarithms [|Tg]; i.e. one should expand in 
powers of l/Nc and not in powers of In (M/m), with M ^ m two widely separated 
scales. Large logarithms are in fact the main source of complications in low-energy flavour- 
changing processes, because the electroweak scale Mw where the short-distance quark 
transition takes place is much larger than the long-distance hadronic scale. 

The large short-distance logarithms can be summed up with the use of the Operator 



Product Expansion (OPE) [|T9| and the renormalization group [^. The proper way to 



proceed makes use of modern Effective Field Theory techniques ||2T|]. One starts above 
the electroweak scale where the flavour-changing process, in terms of quarks, leptons and 
gauge bosons, can be analyzed within the usual gauge-coupling perturbative expansion in 
a rather straightforward way. The renormalization group is used to evolve down in energy 
from the scale Mz, where the top quark and the Z and bosons are integrated out. 
That means that one changes to a different Effective Theory where those heavy particles 
are no longer explicit degrees of freedom. The new Lagrangian contains a tower of operators 
constructed with the light fields only, which scale as powers of 1/Mz- The information 
on the heavy fields is hidden in the (Wilson) coefficients of those operators, which are 
fixed by "matching" the high- and low-energy theories at the point n = Mz- One follows 
the evolution further to lower energies, using the Effective Theory renormalization group 
equations, until a new particle threshold is encountered. Then, the whole procedure of 
integrating the new heavy scale and matching to another Effective Field Theory starts 
again. In this way, one proceeds down to scales fi < rric. 

In this picture, the physics is described by a chain of different Effective Field Theories, 
with different particle content, which match each other at the corresponding boundary 
(heavy threshold). This procedure permits to perform an explicit summation of large 
logarithms t = In (M/m), where M and m refer to any scales appearing in the evolution. 



One gets finally an effective AS = 1 Lagrangian, defined in the three-flavour theory p^^5 



^ei=' = -^K. E C.ii^) Q.(/i) , (4) 
which is a sum of local four-fermion operators Qi, constructed with the light degrees of 
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freedom, 



Q2 = (sm)v_a M)v-A ' 

Q7,9 = l{sd\^^Y.qeq {qq)Y±A ^ 

modulated by Wilson coefficients Ci{fi) which are functions of the heavy masses. Here a, 
P denote colour indices and are the quark charges (e^ = 2/3, = Cs = —1/3). Colour 
indices for the colour singlet operators are omitted. The labels {V ± A) refer to the Dirac 
structures 7^(1 ± 75). 

We have explicitly factored out the Fermi coupling Gp and the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements Vij containing the usual Cabibbo suppression of K de- 
cays. The unitarity of the CKM matrix allows to write the Wilson coefficients in the 
form 

difi) = Zi{fi) + r Viifi) , (6) 

where r = —V^^V^l/V^^V*^. The CP-violating decay amplitudes are proportional to the i/i 
components. 

The overall renormalization scale n separates the short- (M > /i) and long- (m < /i) 
distance contributions, which are contained in Ci{fi) and Qi, respectively. The physical 
amplitudes are independent of /x; thus, the explicit scale (and scheme) dependence of the 
Wilson coefficients should cancel exactly with the corresponding dependence of the Qi 
matrix elements between on-shell states. 

Our knowledge of AS = 1 transitions has improved qualitatively in recent years, thanks 
to the completion of the next-to-leading logarithmic-order calculation of the Wilson coeffi- 
cients [^,|2^. All gluonic corrections of 0{a^t^) and 0{a'^~^^t^) are already known. More- 
over the full rrit/Mw dependence (to ffist order in and a) has been taken into account 
at the electroweak scale. We will fully use this information up to scales /i ~ (9(1 GeV), 
without making any unnecessary expansion in powers of I/Nq- 

In order to predict physical amplitudes one is still confronted with the calculation of 
hadronic matrix elements of quark operators. This is a very difficult problem, which so far 
remains unsolved. As indicated in figure [l|, below the resonance region one can use global 
symmetry considerations to define another Effective Field Theory in terms of the QCD 
Goldstone bosons (vr, K, if). The Chiral Perturbation Theory (xPT) formulation of the 
Standard Model [P^-P^ is an ideal framework to describe the pseudoscalar-octet dynamics, 
through a perturbative expansion in powers of momenta and light quark masses over the 
chiral symmetry breaking scale (A-^ ~ 1 GeV). Chiral symmetry fixes the allowed xPT 
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Figure 1: Evolution from Mw to Mk- 



operators, at a given order in momenta. The only remaining problem is then the calculation 
of the corresponding chiral couplings from the effective short-distance Lagrangian; this 
requires to perform the matching between the two Effective Field Theories. 

It is here where the I/Nq expansion proves to be useful. At leading order in I/Nq, 
the matching between the 3-flavour quark theory and xPT can be done exactly. We will 
determine the needed chiral couplings in the large- A^^c limit, in a quite straightforward way. 
The scale and scheme dependences of the short-distance Wilson coefficients are of course 
completely removed in the matching process, at leading order in 1/Nc- Any remaining 
dependences are higher-order in the 1/Nc expansion and, thus, numerically suppressed; 
they are included in our estimated theoretical uncertainty. 

There is still an important source of large logarithms that needs to be identified and kept 
under control. The FSI of the pseudo-Goldstone pions generate large infrared logarithms, 
involving the light pion mass, which are next-to- leading in 1/Nc- These chiral logarithms 
can be computed within the effective xPT framework. Moreover, as shown in refs. 
they can be exponentiated to all orders in the momentum expansion. Since this is a 1/Nc 
suppressed (but numerically large) effect, it generates an important correction, not included 
in the previous leading-order determination of chiral couplings. 

The paper is organized as follows. The usual isospin formalism for K ^ tttt decays 
and the relevant formulae for e'/ e are collected in section 0. Section ^ presents the low- 
energy xPT description. The matching between the short- and long-distance effective 
theories is performed in section ^, at leading order in 1/Nc- Section ^ summarizes the 
large-A"(7 predictions for the different isospin amplitudes. The one-loop chiral corrections 
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are discussed in section |^. Section |^ incorporates higher-order corrections induced by FSI, 
within the chiral framework. The Standard Model prediction for e' je is worked out in 
section ^ where two different numerical analyses are presented. The first one incorporates 
the experimental value of e, while in the second one its theoretical prediction is used instead. 
Both analyses give compatible results. Our conclusions are finally given in section ^. 
We have collected in several appendices the analytical results from the one-loop chiral 
calculation of the different K tttc amplitudes. 



2 K ^ nn Amplitudes 

We adopt the usual isospin decomposition: 



TT TT 



TT TT 



1 . 

A0-V2A2. 



(7) 



The complete amplitudes Aj = Aj expji^Qj include the strong phase shifts 5q. The S- 
wave TT-vr scattering generates a large phase-shift difference between the 1 = and 1 = 2 



partial waves 33 



(5° - 61) (Mi) = 45° ± 6° 



(8) 



There is a corresponding dispersive FSI effect in the moduli of the isospin amplitudes, 
because the real and imaginary parts are related by analyticity and unitarity. The presence 
of such a large phase-shift difference clearly signals an important FSI contribution to Aj. 
In terms of the K ^ tth isospin amplitudes, 



— = e 
e 



1$ 



Im(A2) Im(Ao) 



Re{A2) Re(Ao) 



Owing to the well-known "A/ =1/2 rule", e'/e is suppressed by the ratio 

u = Re(v42)/Re(Ao) ^ 1/22. 
The phases of e' and e turn out to be nearly equal: 



^2 - AO 



TT 

4 



0. 



(9) 



(10) 
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The CP-conserving amplitudes Re(y4/), their ratio uj and \e\ are usually set to their ex- 
perimentally determined values. A theoretical calculation is then only needed for lm{Aj). 

Using the short-distance Lagrangian (^, the CP-violating ratio e'/e can be written 
as i 



UJ 



--Im(TOe--^ 



IB 



1 p{2) 



00 



[12) 
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where the quantities 

^^'^ = E (M/IQ.I^) (13) 

i 

contain the contributions from hadronic matrix elements with isospin / and 

= (14) 

parameterizes isospin breaking corrections. The factor 1/uj enhances the relative weight of 
the 1 = 2 contributions. 

The hadronic matrix elements {{7nT)i\Qi\K) are usually parameterized in terms of the 
so-called bag parameters Bi, which measure them in units of their vacuum insertion ap- 
proximation values. In the Standard Model, P^^^ and P'-^-* turn out to be dominated by the 
contributions from the QCD penguin operator Qq and the electroweak penguin operator 
Qs, respectively 0]. Thus, to a very good approximation, e' /e can be written (up to global 
factors) as 

s' r„n/2i/. ^ X „ . „C3/2~n ^^^^ 



e 

The isospin-breaking correction coming from vr^-r/ mixing was originally estimated to 
be fl"}^ = 0.25 [^^. Together with the usual ansatz ~ 1, this produces a large 



numerical cancellation in eq. ([T5| ) leading to low values of e'/e around 7-10 ^. A recent 
improved calculation of TT^-rj mixing at 0{p'^) in xPT has found the result p6| 



= 0.16 ±0.03. (16) 
This smaller number, slightly increases the naive estimate of e'/e. 



3 Chiral Perturbation Theory Description 



In the limit rriu, riid, — > 0, the QCD Lagrangian for light quarks has a SU(3)l<^ SU(3)ji 
symmetry, which is spontaneously broken to SU{3)v- The lightest particles of the hadronic 
spectrum, the pseudoscalar octet (tt, K, r/), can be identified with the corresponding 
Goldstone bosons. Their low-energy interactions can be analyzed within xPT p8|^2|, 
which is an expansion in terms of momenta and meson (quark) masses. The Goldstone 
fields are parameterized as 



V 



2" 



in 



TT 

K- 



TT 



2 



\ 



(17) 



h J 



and appear in the Lagrangian via the exponential representation U = exp{\/2i^ / f) , with 
f = 92.4 MeV the pion decay constant at lowest order. Under a chiral transformation 
g = (ql, gji) e SU{3)l ® SU{3)r, the matrix U changes as U gnU g\. 
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The effect of strangeness-changing non-leptonic weak interactions with AS" = 1 is 
incorporated |^ in the low-energy chiral theory as a perturbation to the strong effective 



Lagrangian. At lowest order, the most general effective bosonic Lagrangian, with the same 
SU{3)l (S> SU{3)ji transformation properties and quantum numbers as the short-distance 
Lagrangian (Hf), contains three terms: 



^2'=' = -^K.Cr|^?8 [{XL,L^)+e'f%^{XU^QU) 



+ 927 (^L^23Lii + -L^21-^^13^ I 



where the matrix = —iU'^D^U represents the octet of V — A currents, at lowest order 
in derivatives, Q = diag(|, — |, — |) is the quark charge matrix, A = (A^ — iX'^)/2 projects 
onto the s ^ d transition [Ajj = ^is^j^] and (A) denotes the flavour trace of A. 

The chiral couplings gg and (727 measure the strength of the two parts of the effective 
Lagrangian (|) transforming as (8^,, l/j) and (27^, respectively, under chiral rotations. 
Chiral symmetry forces the lowest-order Lagrangian to contain at least two derivatives 
(Goldstone bosons are free particles at zero momenta). In the presence of electroweak 
interactions, however, the explicit breaking of chiral symmetry generated by the quark 
charge matrix Q induces the C(p°) operator {XWQU) [^^, transforming as (8l,8r) 



under the chiral group. In the usual chiral counting ~ 0{p^) and, therefore, the Qew 
term appears at the same order in the derivative expansion than (^g and (727- One additional 
term |^D| proportional to the quark mass matrix, which transforms as (8^,1/2), has not 
been written in the lowest-order Lagrangian (p!8[) , since it does not contribute^ to physical 
K ^ TTTT matrix elements 

The tree-level K ^ titt amplitudes generated by the 0{p'^) xPT Lagrangian (^) are: 

^0 = KdC v^/ { (^?8 + ^ 927) {Ml - Ml) -p'e' gs ^?e.} , 

^2 = KdK*. If {5 927 {Ml -Ml) -3 f gs g^u.] . (19) 

The strong phase shifts are zero at lowest order. Taking the measured phase shifts into 
account, the moduli of gs and (727 can be extracted from the CP-conserving K ^ 2tt decay 
rates. A lowest-order phenomenological analysis neglecting the tiny electroweak 
corrections proportional to e^gew, gives: 

1^781^5.1, 1^7271^ 0.29. (20) 

The huge difference between these two couplings shows the well-known enhancement of 
octet |AJ| = 1/2 transitions. 

^ The contributions of this term to K ^ tttt amphtudes vanish at 0{p'^), while at 0(p'^) they can be 
reabsorbed through a redefinition of the local 0{p^) AS* = 1 chiral couplings p]-^. 

^ A general analysis of isospin breaking and electromagnetic corrections to ^ tttt transitions is 
presented in refs. p^-|4^. 
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The isospin amplitudes Ai have been computed up to next-to-leading order in the 
chiral expansion [^,^,|45|-^. Decomposing the isoscalar amplitudes in their octet and 
27-plet components as = -^o + -^o 5 the results of those calculations can be written 
in the form: 



(8) 



fl [g^^ (1 + + Ac4 



{ew) 



(21) 



for the octet isoscalar amplitude, 



.(27) 



-fu^VL ^ /. 9.7 {Ml - Ml) [1 + A,4"^ + Ac4"^ 



for the 27-plet isoscalar amplitude and 



f22) 



A = -^K^KTs \ fn {5 927 {Ml - Ml) [1 + AlAT^ + AcAT^ 



- 3 figs [ge^ (1 + Ai4^")) + 



(eui) 
2 



(23) 



for the 1 = 2 amplitude. The electroweak penguin contributions have been also included. 
These formulae contain chiral loop corrections A^^^"* , coming from the lowest-order La- 
grangian ([T8|) and its strong counterpart. Loop corrections are always subleading in the 
1/Nc expansion, so that they do not enter the large- Ai'c matching procedure outlined in 
the introduction. One-loop corrections to K ^ nir have been extensively analyzed in 
ref. [0, with the aim of identifying and resum FSI effects. Those effects, subleading in 
l/Nc but numerically relevant, will be taken into account in sections ^ and |^. 

At next-to-leading order in the chiral expansion, i.e. 0{Gpp'^) and C(Gire^j9^), the 
complete Lagrangian which mediates non-leptonic weak interactions with AS* = 1 can be 
written as follows ||,||,|5|-||]: 



rAs=i 



f (gs E A. of + g27 E A of + gs e^f E or 



(24) 



For the octet and 27-plet weak operators Of and Of the basis constructed in ref. has 
been adopted]^. For the electroweak operators Of^ we use the basisQ of ref. [^]. We refer 
to those references for the explicit form of the operators. 

^ For the octet operators one can use either the basis of ref. or the basis of ref. [Q . For completeness 
we provide the transformation rules between the two bases in appendix 

■* Our operators Of ^ are denoted with Qi in ref. Q and their coupling G% is related to our via the 
identity Gg - - {GfI^) VudVL 9s- 
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The O(p^) and (9(e^p^) tree-level contributions to the K ^ titi amplitudes are easily 
computed with the Lagrangian and its strong counterpart. The complete expressions 
can also be obtained from refs. |Q and [|5|: 



2 A/f ^ 

fn 

+ —f i-2E, - 4^2 - 2^3 + 2E,o + Eu + ^E.^) , (25) 
Ac4''^ = Ac + ^ (D.-D,- 9De + ADj) 

J IT 

+ -jf (-6^1 - 2^2 + 2/^4 + 6^6 + Dr) , (26) 

Ac4''"^ = ^e.Ag-) + ^(Zi + 2Z2) + ^(4Zi + 2^2-^6) (27) 
Ml- - 

- ^.2 " (8^3 - 24^4 + 9^5 + 6Zr - SZg - SZg - 2Zio + 2Zn + 2Z12) , 

Ac4''^ = Ac + ^ iD,-D, + ADj) + ^ (-2D2 + 2D4 + Dj) , (28) 

AcAt^ = ge^ Ag'") + ^ (2Zi + 4^2 - ^e) + ^ (2Zi + Z,) 
Ml - 

+ 0^9 " (-4^3 + 12^4 - 3^8 - 3Z9 - 2Zio + 2Zu + 2Z12) , (29) 



3/^ 



where 



Ac = -^(Mi + 3M^^)-1^(2M^ + M^ 

J TT J TT 

= {Ml + 5M^^) - ^ (2Mi + M^) . (30) 

J 7T J TV 



There are seven (8^,, 1/j) operators Of {i = 1, 2, 3, 10, 11, 13, 15), six (27^, 1r) operators 



Of {i = 1,2,4,5,6,7) and twelve electroweak operators Of (z = 1, . . . ,12) contribut- 
ing to — > TTTT matrix elements P3|,^. The practical limitation of a systematic xPT 
evaluation of the K ^ -k-k isospin amplitudes is in the fact that the counterterms which 
appear at next-to-leading order are not fully known and their determination would require 
the experimental knowledge of a large set of weak AS = 1 processes. 

In addition, there are contributions involving the lowest-order AS" = 1 Lagrangian (|IH]) 
combined with the 0{p'^) strong chiral operators with couplings L,, introduced in ref. p9[ . 
In previous analyses [^^^0 these corrections, shown in eqs. (PPD, were factorized 
as global factors in front of the corresponding amplitudes: 1 + Aq = f'^ / {flfx) ~ 0.65, 
1 + A^^"''' = /{f^fx) ~ 0.58. A factor f^/{f^fK) arises from wave-function renormali- 
zation, while the remaining powers of f/fn are needed to rewrite in terms of the physical 
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pion decay constant the explicit dependences of the tree-level amplitudes (|1^) on the chiral 
Lagrangian coupling /. This procedure induces a sizeable suppression which is finally 
compensated by large and positive corrections from the 0{p^) weak counterterms. We 
prefer to keep all local contributions together and perform a consistent large-A^^c 

calculation of their global size. 



4 harge—Nc Matching 

In the large-Nc limit the T-product of two colour-singlet quark currents factorizes: 

(J -J) = (J) (J) + . (31) 

In other words, colour exchanges between the two currents J are I/Nq suppressed and 
in this limit the factorization of four-quark operators is exact. Since quark currents have 
well-known realizations in xPT |29H3l|, the hadronization of the weak operators Qi can 
then be done in a quite straightforward way. Thus, at large-A^c the matching between 
the short-distance Lagrangian (1) and its long-distance xPT realization can be explicitly 
performed. 

The chiral couplings of the lowest-order Lagrangian (|18D have the following large-iVc 
values: 

gT = -^gi(/^) + lcM + c,{f,) - 16L5 (^ ^^^^^^ ^ CM, 
gTr = llcM + CM], (32) 
igse^g..r = -3(M^)^3(^). 

Together with the C(p^) amplitudes in eqs. (|19]), these results are equivalent to the 



standard large- A^^c evaluation of the usual bag parameters Bi. In particular, for e'/e, 
where only the imaginary part of the Qi couplings matter [i.e. Im(Cj)], eqs. (^) amount 
to B^^'^'^ ^ B^^'^^ = 1. Therefore, up to minor variations on some input parameters, 
the corresponding e'/e prediction, obtained at lowest order in both the 1/Nc and xPT 
expansions, reproduces the published results of the Munich and Rome p groups. 

The large-iVc limit has been only applied to the matching between the 3-fiavour quark 
theory and xPT, as indicated in figure |I]. The evolution from the electroweak scale down to 
fi < rric has to be done without any unnecessary expansion in powers of I/Nq] otherwise, 
one would miss large corrections of the form In (M/ m) , with M ^ m two widely 
separated scales [0. Thus, the Wilson coefficients contain the full dependence. 

The operators Qi {i 7^ 6,8) factorize into products of left- and right-handed vector 
currents, which are renormalization-invariant quantities. The matrix element of each single 
current represents a physical observable which can be directly measured; its xPT realization 
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just provides a low-energy expansion in powers of masses and momenta. Thus, the large- 
Nc factorization of these operators does not generate any scale dependence. Since the 
anomalous dimensions of Qi {i 7^ 6,8) vanish when Nc —>■ 00 [1^, a very important 
ingredient is lost in this limit [Q. To achieve a reliable expansion in powers of 1/Nc, one 
needs to go to the next order where this physics is captured l5y,|51[ . This is the reason why 
the study of the AI = 1/2 rule has proved to be so difficult. Fortunately, these operators 
are numerically suppressed in the e' /e prediction. 

The only anomalous dimensions which survive when Nq — »• 00 are the ones corre- 
sponding to Qq and Qs ll^,|3^. One can then expect that the matrix elements of these 
two operators are well approximated by this limitQ [pO|-|52[|. These operators factorize into 
colour-singlet scalar and pseudoscalar currents, which are dependent. Since the products 
nig g(l, 75)g, are physical observables, the scalar and pseudoscalar currents depend on /i 
like the inverse of a quark mass. Conversely, the Wilson coefficients of the operators Qq 
and Qs scale with /i like the square of a quark mass in the large- A'^p limit. 

The xPT evaluation of the scalar and pseudoscalar currents provides, of course, the 
right /i dependence, since only physical observables can be realized in the low-energy 
theory. What one actually finds is the chiral realization of the renormalization-invariant 
products nig 5(1,75)^. This generates the factors 



nig = 



P 



Bp _ _-Bo /tt 

' / " U f 

Ml 

1 + 



J n 



2Ml + Ml, , Ml , 

72 ^(3^4 - 4L6) - 8 -f (2L8 - L5) 



3 I/^ - Z/K - - Z^r, 



(m^ + nig) (/i) 



1 + 



4L5 
P 

J n 

M 



Ml + 4 



2Ml + Ml 

J TT 



(3L4 - 4L6) 



K 



P 

J TT 



(2^8 - L5) - 2 z/^ - z/x - - z/^ 



(33) 



in eqs. (^21), which exactly cancel the dependence of C^^siP) at large Nc [|T8i^ ,^^-|52|. It 
remains a dependence at next-to-leading order. The parameter Bq is a low-energy coupling 
of the 0{p'^) strong chiral Lagrangian, which accounts for the vacuum quark condensate 
at lowest order in the momentum expansion. The one-loop corrections vp {P = it, K, 77), 
defined in appendix are identically zero in the limit Nc 00. 

While the real part of gs gets its main contribution from C2, lm{gs) and lm{gsgew) 
are governed by Cq and Cs, respectively. Thus, the analyses of the CP-conserving and 
CP-violating amplitudes are very different. There are large 1/Nc corrections to R.e{gi) 

^ Some insight on these matrix elements can be obtained from the two-point fmictions '^ii{q^) = 
i J d^x e**^ {T{Qi{x)Qi{Oy)) , since their absorptive parts correspond to an inclusive sum of hadronic matrix 
elements squared. The known 0{as) results [^-p2| show that the large- A^c limit provides an excellent 
approximation to 5*66: but an incorrect description of ^'22- 
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)2|, which are needed to understand the observed enhancement of the (8l, 1r) couphng. 
However, the laige-Nc hmit can be expected to give a good estimate of lm{gi). 

Contrary to the other Qi operators, the leading-order contribution of Qq involves the 
coupling L5 of the 0{p'^) strong chiral Lagrangian. The large-Xc value of this chiral 
coupling can be estimated from the ratio of the kaon and pion decay constants: 



T 00 

-^5 




2.1 • 10" 



(34) 



The Qq contribution dominates the numerical value of lm{g^). In the large-Nc limit, the 
combined effect of all other operators only amounts to a 5% correction. 

The 0{p^) corrections introduce dependences on three additional strong chiral cou- 
plings. At large Nc, 

= = 0. (35) 

To determine Lg, we impose the stronger requirement of lowest-meson dominance [^ ,|5^ 
and assume that the scalar form factors vanish at infinite momentum transfer. This implies 



the relation 55 



(2Lg 







(36) 



which is well satisfied by the phenomenological determinations of those constants [P9|,|5C 

The operators Q3 and Q5 start to contribute at O(p^), while the electroweak penguin 
operators Q7, Qg and Qio give their first contributions at (9(e^p^). The large-iVc matching 
at the next-to-leading chiral order fixes the couplings Ei, Di and Zi of the long-distance 



chiral Lagrangian (p^). We only quote the values of those couplings contributing to K 
717T amplitudes. 

For the 0{p'^) couplings, one gets: 



{gsEsY 
(gsEioY 



{gsE 



-48X,,(M^) CM 




2 L, 



5 5 



—8 (2X14 + 2X15 + X38 
4L5 

-16 (Xi5 + 2X17 - X38 



2 3 

5 5 




CM 



12 



13; 



(^8^15)^ 



:-X3. + X3s) {^^^\cM 



P 



(37) 



All other {27 1^) couplings contributing to i^' tttt -D5, -De and 1^7) are zero 

at large- Ai'c. 

The 0{p'^) contributions from the operator have been computed using the 0{p^) 
Lagrangian of ref. the couplings Xi refer to the list of 0{p^) SU{3) operators given in 
that reference. These couplings however are unknown, so in practice the Qq contribution 
is missing in eqs. (^7\). The remaining terms are in agreement with the results obtained in 



ref. [48 



The non-zero 0{e p ) couplings relevant for K ^ nn are: 





^ -24 




{gse'z.r 
[gse'z.r 


= C'io(^), 




( 2 ry \00 
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5 Isospin Amplitudes at Leading Order in l/N( 



c 



Combining the results of the previous sections, one gets the predicted K amplitudes 
at leading order in 1/Nc- The different contributions to the isospin amplitudes take the 
following form: 



1 + Ac4'^l" 

+ ^ C2(/x) + C4(/x) - 16 L5 CM 



K 



7r) 5 



(39) 



1 + = g^, 



1 + Ac4'^)l°" 



[Ciif^) + C2if^)]ftiM-^ 



(40) 



13 



-3C8(/i) 



Mi 



{rus + mg)(/i) 



1 + -:#M2 



J 7r 



- [Cj - Cg 



Ml -Ml 



CM ^-^^^^ ft {Ml) 

J TT 



(41) 



-3C8(/i) 



Mi 



1 + ^Ml 



P 

J n 



+ 2[^7-C9-Cio](/i) 



M2 



M 



J TT 



(42) 



For the operators Qi {i 7^ 6,8), which are products of colour-singlet vector and axial- 
vector currents, these are exact large- A^^c results to all orders in the chiral expansion, as can 
be easily seen factorizing the operators at the quark level. The xPT framework discussed 
before reproduces these results in a perturbative way, through the momentum expansion 
of the Kn scalar form factor at 0{p^): 



ft {Ml) ^ ft {Ml) + 



Ml 



f^-{Ml) 



Mi - Ml 



1 + 



4L5 n 
J IT 



(43) 



The form factors f±^{t) are defined through the matrix element of the vector current, 



(7r|s7'^g|i^) = Ck. {{Pk + K)" /^'^(t) + (P/^ - M'^ /"""W} 



{q = u,d), (44) 



where t = {Pk — ^tt)^, CxOnO = — Cx+ttO = and CxOn- = Ck+-k+ = — 1- 

The wave-function renormalization corrections Ac [eq. (PPj)] have been cancelled by 
weak 0{p^) contributions, as it should since we are dealing with conserved currents. Once 
the 0{p'^) results are written in terms of the physical pion decay constant /tt, higher-order 
chiral contributions only introduce the small correction factor ft {Ml) ^ 1.02. 

The hadronic matrix elements of the operators Qq and Qs factorize into products of 
scalar and pseudoscalar currents, which cannot be directly measured. The xPT predictions 
are then needed to determine those hadronic currents. The electroweak penguin matrix 
elements are known to O(p^). Again, one observes that the contributions from local weak 
terms {Zi and Zq) cancel the negative contribution from A^^^ and reverse the sign of the 
0{p^) correction. The contribution of the penguin operator Qq is only known at 0{p'^). 
For Qq we cannot just include the Ac correction, because the corresponding weak O(p^) 
counterterms are unknown and large cancellations can be expected. In eq. (|39|) we have 
taken a global correction factor f^'^{Ml) for the octet amplitude. This is a reasonable 
assumption^ since nearly all known pieces have this common correction. Only {Qs)o gets 
a different (and larger) correction. 

^ In fact, the factor f^'"{Ml) already appears in the lowest-order Qq contribution to gg, through the 
0{p'^) correction in eq. (|3|). 
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Table 1: Numerical values of the weak cliiral couplings in the large- Ai'c limit. 





"1 + AcA^^^' 


00 




(1.3 ±0.2 ±0.1) + r (1.12 ± 0.08 IgJ^) 




I + Ac At' 


CO 




(0.47 ±0.01 ±0.0) 




'g,^ + AcAi'-y 


00 


- (0.085 ± 0.085 1[J:[}|) - r (2.33 ± 0.07 


e'gr 


''ge^ + AcAt^'' 


00 


- (0.07 ± 0.07i[J;|]^) - T (1.34 ± 0.03 toit) 



The scalar and pseudoscalar currents introduce a quadratic dependence on quark masses 
in the contributions from the operators Qq and Qs- At present, the most reliable de- 
terminations of the light quark masses give ms{l GeV) = (150 ± 25) MeV [p^-|63| and 
(m„ ± md){l GeV) = (12.8 ± 2.5) MeV 0, at the scale = 1 GeV. We then take: 



m, 



±mg)(l GeV) = (156 ±25) MeV. 



(45) 



Table shows the resulting numerical predictions for the weak chiral couplings. The 
central values have been obtained at = 1 GeV. The first errors indicate the sensitivity 
to changes of the short-distance renormalization scale in the range Mp < /i < rric and to 
the choice of 75 scheme in the next-to-leading order calculation of the Wilson coefficients. 
The second uncertainties correspond to the input values of the quark masses. 

For historical reasons, the values of the short-distance Wilson coefficients are usually 
given in terms of Aqcd (in the three or four flavour theory). Nowadays, that as is experi- 
mentally known with rather good accuracy, it is unnecessary to introduce this additional 
auxiliary parameter which only complicates the flnal expressions. Since the most impor- 
tant as corrections appear at the lowest scale /i ~ 0{1 GeV), we have flxed the strong 
coupling at the r mass, where it is known with about a few percent level of accuracy: 



as{Mr) = 0.345 ±0.020 . 



(46) 



The high-energy matching scale is chosen to be intermediate between the IV-boson and 
the top quark mass scale. We have performed the matching directly at the Z'-boson mass 



scale where ag is best known [|66| , 



as{Mz 



0.119 ±0.002 



(47) 



The measured values (^) and (^) are in perfect agreement, if one performs ||6^ a four- 
loop evolution of as between Mz and Mr, with the appropriate matching conditions at 
the different thresholds [^]. The values of as at the other needed scales can be deduced 
from (1^). The numerical uncertainties associated with the present error on as have been 
included in our results, but they are negligible in comparison with the uncertainties from 
other sources. 

The dominance of Qq and Qs in the CP-odd amplitudes (the ones proportional to the 
CKM factor r) is apparent in table where those pieces show a very strong dependence 
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on quark masses (second error bars). In comparison, the short-distance uncertainties are 
much smaller. The opposite behaviour is observed in the CP-conserving couplings Re^gs) 
and Re(5'27), which are dominated by Qi and Q2. The 27-plet coupling, which does not 
get any penguin contribution, satisfies lm{g27) = for all practical purposes. 

Taking fl^^ = 0.16, Im(V^*V^^) = 1.2 ■ 10~^ and the central values in table |l] for the 
CP-odd amplitudes, one gets the \a.Tge-Nc prediction Re{e'/e) ~ 0.8 ■ 10^'^. Although 
numerically suppressed, the operators Qi, Q2 and Q4, which are not well approximated by 
the large- A'^c limit, provide also small corrections to Im(y4o). In refs. the measured 

CP-conserving rates are used to estimate those contributions. This amounts to multiply 
the corrections from these operators by a factor ^0 ~ 4.9, to compensate for the under- 
estimated coupling Re{gs). Adopting this prescription, one gets Re{e'/e) ~ 0.5 ■ 10~^, in 
agreement with the findings of refs. 



6 Chiral Loop Corrections 

The previous tree-level amplitudes do not contain any strong phases Sq. Those phases 
originate in the final rescattering of the two pions and, therefore, are generated by chiral 
loops which are of higher order in the 1/A"c expansion. Analyticity and unitarity require the 
presence of a corresponding dispersive FSI effect in the moduli of the isospin amplitudes. 
Since the strong phases are quite large, specially in the isospin-zero case, one should expect 
large higher-order unitarity corrections. Intuitively, the behaviour of the / = and 1 = 2 
S-wave phase shifts as a function of the total invariant mass of the two pions suggests a 
large enhancement of the 1 = amplitude and a small suppression of the 1 = 2 amplitude. 

The size of the FSI effect can be estimated at one loop in xPT. The dominant one-loop 
correction to the octet amplitude comes indeed from the elastic soft rescattering of the 



two pions in the final state. The existing one-loop analyses |42, 53| show that pion loop 
diagrams provide an important enhancement of the Aq amplitude by about 40%, implying 
a sizeable reduction of the phenomenologically fitted value of Igg] in eq. ([20|). 



The complete formulae for the one- loop corrections A^A) are compiled in the ap- 
pendices. The usual one-loop function _B(Mf,M|,p^) is defined in appendix while 
appendix y contains explicit results for the different isospin amplitudes. The contribu- 
tions proportional to B{Mj,, Mj>, M]^), with P = ir, K, r], arise from intermediate vrvr, 
KK and rjr] states. At s = (p^^ + = ^^e only possible absorptive contribution 
comes from the elastic vrvr rescattering: 

Ai4^) = -i (2M1. - M^) B{Ml Ml Ml) + . . . (48) 



Al4''^ = \(mI-2MI) B{MlMlMl) + ... (49) 



2 
1 

2 

where 

B( Ml Ml Ml) = —1— 
V KJ (47r/^)2 1 " 



1 / 1 + ^vr 1 

m I \ — tir 

1 — a. 




(50) 
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with u the chiral loop renormahzation scale and 



= Jl-—f. (51) 



Thus, all isoscalar amplitudes get the same absorptive contribution, as it should, since 
they have identical strong phase shifts. The same is true for the two amplitudes with 
1 = 2. The one- loop absorptive contributions reproduce the leading xPT values of the 
strong rescattering phases 6q, with / = 0, 2: 

tan S','\Ml) = (2MI - Ml ; 2MI - M^) . (52) 

The numerical values of 5o'^(M|;) predicted by xPT at leading order, 5q[M]^) = 25° and 
6q{M'^) = —12°, are significantly lower than their experimental values, implying that 
higher-order rescattering contributions are numerically relevant. The phase-shift differ- 
ence, ^0 — ^0 = 37°, is slightly less sensitive to higher-order chiral corrections . 

The 2-71 intermediate state induces a large one-loop correction to the 1 = amplitudes. 
At u = Mp, the 27r contribution to the isoscalar amplitudes is A^^q = 0.43 + 0.46 i, 
while AiA2^''\mT = —(0.19 + 0.20 z); i.e. the expected enhancement (suppression) of the 
J = (/ = 2) amplitudes. The contributions from other one-loop diagrams, not related to 
FSI, are different for the different amplitudes Aj . 

Let us write our isospin amplitudes in the form 



where 4^^°° the large-A'^c results obtained in the previous section. The correction 
factors C) contain the chiral loop contributions we are interested in. At the one-loop 
level, they take the following numerical values: 

C^Q^ ^ I + AlA^ = 1.27 ±0.05 + 0.46 z, 

Cj^^^ ^ I + AlA'^'^^ = 2.0 ±0.7 + 0.46 z, 

Cf""^ ^ 1 + AlX^"'^ = 1-27 ±0.05 + 0.46 z, 

^ 1 + Al4^^^ = 0.96 ± 0.05 - 0.20 i , 

4'"^ ^ 1 + Ai4'"^ = 0.50 + 0.24- 0.20 ^. (54) 

The central values have been evaluated at the chiral renormahzation scale u = Mp. To 
estimate the corresponding uncertainties we have allowed the scale u to change between 
0.6 and 1 GeV. The scale dependence is only present in the dispersive contributions and 
should cancel with the corresponding u dependence of the local counterterms. However, 
this dependence is next-to- leading in I/Nq and, therefore, is not included in our large- A^^^ 
estimate of the 0{p^) and (9(e^p^) chiral couplings. The u dependence of the chiral loops 
would be cancelled by the unknown l/A'c'-suppressed corrections AcA^j^\i') — AcA\^^°°, 
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that we are neglecting in the factors C) . The numerical sensitivity of our results to the 
scale z/ gives then a good estimate of those missing contributions. 

The absorptive contribution induces a large one-loop correction to the / = amplitudes. 

(07) 

The dispersive correction to A^^q is even larger, but it has a smaller phenomenological 
impact because the isoscalar K nn amplitude is dominated by its octet component; this 
27-plet correction has a strong dependence on u and, therefore, a rather large uncertainty. 
Although the one-loop correction to the 1 = 2 (27^, Ir) amplitude is rather moderate, the 
electroweak 1 = 2 amplitude gets a large dispersive correction with negative sign. This 
induces a corresponding suppression of l^'""^! by about 46%. 

The numerical corrections to the 27-plet amplitudes do not have much phenomenolo- 
gical interest for CP-violating observables, because Im((727) = 0. Remember that the CP- 
conserving amplitudes Re(A/) are set to their experimentally determined values. What is 
relevant for the e'/e prediction is the 35% enhancement of the isoscalar octet amplitude 
Im[y4o'^''] and the 46% reduction of Im[y42'^'"'']. Just looking to the simplified formula ([l5|), 
one realizes immediately the obvious impact of these one-loop chiral corrections, which 
destroy the accidental lowest-order cancellation between the / = and 1 = 2 contributions, 
generating a sizeable enhancement of e'/e. 

A complete 0{p'^) calculation |^,|^ of the isospin-breaking parameter Qjb is not yet 



available. The value 0.16 quoted in eq. (|TB|) only accounts for the contribution from n^-r] 
mixing |]5B| and should be corrected by the effect of chiral loops. Since \C^'^^\ ~ 0.98 ±0.05, 
one does not expect any large correction of Im(A2)/B, while we know that Im[74o^''] gets 
enhanced by a factor 1.35. Taking this into account, one gets the corrected value 



c. 



(27) 



(8) 



0.12 ±0.05, (55) 



where the quoted error is an educated theoretical guess. This value agrees with the result 
^iB = 0.08 ± 0.05 ± 0.01, obtained in ref. [^] by using three different models []7|,|4^, pl| , |53| . 



|7D|,0 to estimate the relevant 0{p ) chiral couplings. 

The one-loop corrections increase the large-iVc estimate from e'/e ~ 0.8 • 10^^ toQ 
e'/e ~ 1.8 ■ 10"^. The contributions to Im(Ao) from the operators (5i,2,4 can be corrected 
phenomenologically, as advocated in ref. this requires now a smaller factor ~ 3.5, 
which results in0 e'/e ~ 1.5 ■ 10"'^. 



7 Final State Interactions at Higher Orders 

Given the large size of the one-loop contributions, one should worry about higher-order 
chiral corrections. The fact that the one-loop calculation still underestimates the observed 
^0 phase shift indicates that a further enhancement could be expected at higher orders. 

^This number is obtained taking the experimental value for e and Im{Vf*Vf.^) = 1.2-10"''. Using instead 
the theoretical prediction for e, one would get e'/e « 2.2 • 10~^. See section 8 for more details on this 
second kind of numerical analysis. 

^ Using the theoretical value of e, one finds e'/e « 1.8 • 10^'^. 
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The large one-loop FSI correction to the isoscalar amplitudes is generated by large 
infrared chiral logarithms involving the light pion mass [0]. These logarithms are universal, 
i.e. their contribution depends exclusively on the quantum numbers of the two pions in 
the final state 0. As a result, they give the same correction to all isoscalar amplitudes. 
Identical logarithmic contributions appear in the scalar pion form factor ||29[, where they 
completely dominate the O(p^) xPT correction. 

Using analyticity and unitarity constraints , these logarithms can be exponentiated 



to all orders in the chiral expansion The result can be written as: 

Cf ^ Cf\Ml) = nj{Ml,s,)Cf\so). (56) 
The Omnes |72-74| exponential! 



o / \ is'is) <» / N { {s ~ So) f dz 5l{z) 

^i{s,Sq) = e''o('> ^j{s,So) = exp\^ / . 

[ TT J [Z — So) [z — s — te) 



(57) 



provides an evolution of C) (s) from an arbitrary low-energy point Sq to s = (ptti + P 
M|-. The physical amplitudes are of course independent of the subtraction point sq. Intu- 
itively, what the Omnes solution does is to correct a local weak K ^ tttt transition with an 
infinite chain of pion-loop bubbles, incorporating the strong tttt — > tttt rescattering to all 
orders in xPT. The Omnes exponential only sums a particular type of higher-order Feyn- 
man diagrams, related to FSI. Therefore, eq. (^) does not provide the complete result. 
Nevertheless, it allows us to perform a reliable estimate of higher-order effects because it 
does sum the most important corrections.[3 

The Omnes resummation of chiral logarithms is uniquely determined up to a polyno- 
mial (in s) ambiguity, which has been solved with the large- A^^c amplitude Wf^^"" . The 
exponential only sums the elastic rescattering of the final two pions, which is responsible 
for the phase shift. Since the kaon mass is smaller than the inelastic threshold, the virtual 
loop corrections from other intermediate states {K —>■ Kn, Kri,ri7], KK —>■ tttt) can be 
safely estimated at the one loop level; they are included in Cj^\so)- 

Taking the chiral prediction for 6q{z) and expanding Qj{M]^, so) to O(p^), 

(M|-so) f dz 5i{z) 



one should reproduce the one-loop xPT result. This determines the factor C^j^\so) to 
0{p^) in the chiral expansion: 

cf\so) = cfUl-5fi,(M|,so)l ^ l + ALAf-6Qi{Ml,so). (59) 



^ Equivalent expressions with an arbitrary number of subtractions for the dispersive integral can be 
written 

Indeed, a recent dispersive analysis of higher-order corrections induced by Kn rescattering corrobo- 
rates that these additional "crossed-channel" contributions are negligible [75|]. 
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It remains a local ambiguity at higher orders 0,|72|,|7^ . 

Eq. (|56D allows us to improve the one-loop calculation, by taking sq low enough that 
the xPT corrections to Cj^\sq) are moderate and exponentiating the large logarithms 
with the Omnes factor. Moreover, using the experimental phase shifts in the dispersive 
integral one achieves an all-order resummation of FSI effects. The numerical accuracy of 
this exponentiation has been successfully tested through an analysis of the scalar pion 
form factor, which has identical FSI than Aq. 

At So = 0, the dispersive parts of the experimentally determined Omnes exponentials 
are §: 

%(M|;,0) = 1.55 ±0.10, 3?2(M|,0) = 0.92 ±0.03. (60) 

The quoted errors take into account uncertainties in the experimental phase-shifts data 
and additional inelastic contributions above the first inelastic threshold. These numbers 
fit very well with the findings of the chiral one-loop calculation discussed in the previous 
section. The corrections induced by FSI in the moduli of the decay amplitudes Ai generate 
an enhancement of the A J = 1/2 to A J = 3/2 ratio [jl|. 



Ko(M|, 0)/5J2(M|, 0) = 1.68 ± 0.12 



(61) 



This factor multiplies the enhancement already found at short distances. 

At C(p'^), the previous numbers should be corrected with the factors Cj^\so), which 
incorporate additional one-loop contributions not related to FSI. These factors compensate 
the obvious sq dependence of the Omnes exponentials, up to 0{p^) corrections. To estimate 
the remaining sensitivity to this parameter, we have changed the subtraction point between 
So = and So = 3M^ and have included the resulting fluctuations in the final uncertainties. 



The detailed numerical analysis is given in appendix 0. At u 
values for the resummed loop corrections: 



Mp, we get the following 



^(27) 



(ew) 



M27) 

Mew) 
'-2 



3fJo(Mi,so)cS'^(so) 

3fJo(Mi,so)Cr^(so) 
3?2(Mi,so)Cf)(so) 



1.31 ±0.06, 

2.4 ±0.1, 
1.31 ±0.07, 

1.05 ±0.05. 
0.62 ±0.05. 



(62) 



To derive the Omnes representation, one makes use of Time-Reversal invariance, so that 
it can be strictly applied only to CP-conserving amplitudes. Nevertheless, the procedure 
can be directly extended to the CP-violating components relevant for the estimate of e'/ e. 
Working to first order in the Fermi coupling, the CP-odd phase is fully contained in the ra- 
tio of CKM matrix elements r which appears in the short-distance Wilson coefficients and. 



therefore, in J^^^°^ . Decomposing the isospin amplitudes as J^[^^ = AY^''~'^ + t A'j 

the Omnes solution can be derived separately for the two amplitudes A\^^ and A^^'' 
which respect Time-Reversal invariance. 



,{R)CP 
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8 Numerical Analysis 



The CP-violating ratio e' je is proportional to the CKM factor Im(Vj*VJ^). The standard 
unitarity triangle analyses [j70[ have estimated this parameter to be in the range 

Im(V^,:\/J = (1.2±0.2)-10-^ (63) 

This determination is obtained combining the present information on various flavour- 
changing processes; mainly, e, Bq-Bq mixing and the ratio r(6 u)/T{h — > c). The 
final number is sensitive to the input values adopted for several non-perturbative hadronic 
parameters and, thus, there are large theoretical uncertainties which are not easy to 
quantify. 

Since the Standard Electroweak Model has a unique source of CP violation, the same 
combination of CKM factors appears in the theoretical prediction for e, which is propor- 
tional to the K^-K^ matrix element of the AS* = 2 operator: 

m{sLlM{-SLl'dL)\K') ^ ^flMlBK{li) . (64) 

The factor B^ifJ') parameterizes this hadronic matrix element in vacuum insertion units. 
The corresponding Wilson coefficient C/\s=2{^^) is known at the next-to-leading logarithmic 
order [^|7^. Taking appropriate values for the different inputs one finds: 

\e\ = ^BKlm{V,:Vj (18.9 - 14.4 p) , (65) 

with p one of the two CKM parameters, in the Wolfenstein parameterization, which 
characterize the upper vertex of the unitarity triangle. The standard analyses |]76[ favour 
the range p = 0.2 ± 0.1, implying 

\e\ = ^ Bk Im {y:y,,) (16.0 ± 1.4) , (66) 

where 13 ^ = C*as'=2(/^) Bx{p) is the scale-invariant bag parameter. In the large- A^^c limit, 
Bk = BM = 3/4. 

The numerical values of both Im (y^*Vf.j) and p depend on hadronic inputs. However, e 
is rather insensitive to the precise value of p; it changes by less than 10% when p is varied 
within the previously quoted range. 

Thus, we can make two different numerical analyses of e' / e: 



1. The usual one, taking the experimental value of e and adopting the range ( pBD for 
the relevant CKM factor. 

2. Using instead the theoretical prediction of e in eq. (|66|), the ratio e'/e does not depend 
on Im (V^*\^^) p]. The sensitivity of this CKM factor to different hadronic inputs is 
then reduced to the explicit remaining dependence on 13 k- 
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The second type of analysis is more suitable to a systematic 1/Nc approach. The 
theoretical prediction for e' /e depends on ratios of hadronic matrix elements, i.e. Bi/ Ex- 



it is known IjT^ that Bk has sizeable large- A^^c |51,80,81| and chiral |82] corrections, which 
are of opposite sign and could then cancel to some extent. Thus, one can expect the limit 
Nc — > cxo to provide a good starting point to analyze the relevant ratios B^^^"^^ /Bk and 
Bi'/'^/BK. 

We have performed the two types of numerical analysis, obtaining consistent results. 
This allows us to estimate better the theoretical uncertainties, since the two analyses have 
different sensitivity to hadronic inputs. The contributions to Im(Ao) from the operators 
Qi,2,4 have been estimated, following the strategy adopted in ref. i.e. we have corrected 
them with the factor ~ 3.5. 

As a first estimate, we can perform the calculation of e'/e to O(p^) in xPT, without 
making any Omnes resummation of higher-order corrections. Once the large one-loop 
corrections are taken into account, all important ingredients are already caught. We find, 
for the two different types of analysis: 

Reie'/e) = 1.5 ■ IQ-^ ^"^^^f^l^^ = 1.8 ■ 10-^ (67) 

To quantify the uncertainties, we need to consider higher-order effects. Performing the 
Omnes resummation, as indicated in eq. (|56|), one finds: 

Re (e'/e) = lA-W' ^""^^^'^^^^^ = 1.6 -IQ-^. (68) 



These numbers are quite close to the one-loop results (|67|), which indicates that the error 
induced by the chiral loop calculation is not large. 
From the previous numbers, we derive: 

Re (e'/e) = (l.7 ± 0.2 ± 0.5) ■ 10"^ (69) 



The first error indicates the sensitivity to the short-distance renormalization scale, which 
we have taken in the range Mp < fi < rric. The uncertainty coming from varying the 
strange quark mass in the interval (m^ + fng){l GeV) = 156 ± 25 MeV is indicated 



by the second error. We have added a 30% uncertainty from unknown next-to-leading in 
1/Nc local contributions (third error). 



9 Discussion 

The infrared effect of chiral loops generates an important enhancement of the isoscalar 
K ^ irn amplitude. This effect gets amplified in the prediction of e'/e, because at lowest 
order (in both 1/Nc and the chiral expansion) there is an accidental numerical cancellation 
between the / = and 1 = 2 contributions. Since the chiral loop corrections destroy this 
cancellation, the final result for e'/e is dominated by the isoscalar amplitude. Thus, the 
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Standard Model prediction for e' je is finally governed by the matrix element of the gluonic 
penguin operator Qg- 

There are three major ingredients in our theoretical analysis: 

1. A short-distance calculation at the electroweak scale and its renormalization-group 
evolution to the three-flavour theory (/x ^ "nic), which sums the leading ultraviolet 
logarithms. 

2. The matching to the xPT description. 

3. Chiral loop corrections, which induce large infrared logarithms related to FSI. 



The first step is already known at the next-to-leading logarithmic order p6|,p7|. The 
short-distance results are then very reliable. 

We have tried to achieve an acceptable control of the large infrared chiral corrections, 
which are fully known at the one-loop level. A complete two-loop xPT calculation is not 
yet available. Nevertheless, since the leading one-loop corrections are generated by the 
FSI of the two pions, we can use the Omnes resummation to get an idea about the size 
to be expected for the unknown higher-order contributions. The Omnes exponential only 
sums a particular type of higher-order Feynman diagrams, related to FSI. Although it does 
not give the complete result, it allows us to estimate the theoretical uncertainty in a very 
reliable way, because it does sum the most important corrections. The one-loop results 
and the Omnes calculation agree within errors, indicating a good convergence of the chiral 
expansion. 

The most critical step is the matching between the short- and long-distance descrip- 
tions. We have performed this matching at leading order in the 'VjNc expansion, where 
the result is exactly known to Oijy'^ and (9(e^p^) in xPT [O(p^) for Qg]- This can be ex- 
pected to provide a good approximation to the matrix elements of the leading and Qg 
operators. Since all ultraviolet and infrared logarithms have been resummed, our educated 
guess for the theoretical uncertainty associated with l/A^c- corrections is ~ 30%. 

As our final result we quote: 

Re(e7e) = (1.7 ± 0.9) • 10"=^ . (70) 

A better determination of the strange quark mass would allow to reduce the uncertainty 
to the 30% level. In order to get a more accurate prediction, it would be necessary to have 
a good analysis of next-to-leading l/A^c- corrections. This is a very difficult task, but 
progress in this direction can be expected in the next few years 0,|lO|,|l],|8l|,||,|8l. 
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A Octet basis transformation rules 



Following the same notation as the original references, one can change from the octet basis 
J2i Ei Of of ref. |^ to the one of ref. [^], J2i using either the following identities 

for the operators, 



o 



10 ' 



W- 



10 



Of 



wi 



:0 



12 ) 



•-^13 1 



O^ + o^ 



O^ ^O^ , 



08 . 

^15 ) 



W- 
W- 



11 



r)8 . 



12 



W- 



13 



O 



3 1 



4 ) 



w. 



36 



Of-Ol + 



5 1 



(A.l) 



or the coefficient relations 



iV5 


= EiQ — Ell ; 


iVio 


= El — E^ 




= Ell + 2-E'i2 ; 


iVii 


= E2 ; 




= 2-^11 + -^13 ; 


iVl2 


= Es + E, 




= -£^11 ; 


iVl3 


= E4 ; 




= Ei5 ; 




= E,. 



(A.2) 



B One-Loop Functions 

The one-loop function B{Mf,M2,p'^) is defined by the (dimensionally regularized) basic 
scalar integral with two bosonic propagators: 

' I (Cp (Mf-,^) |Mj-(p-,)^i - f' «(^^' • ''"'"'^ + + f^') 

It can be expressed in terms of the function Ji2(p^) p9|l : 



with 



= 16^r-2 ^^Mf + W^^"W 



fB.r 



^^7-^ ^9. 7772 ^72^7 (i^-^) 
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where 

= \\p^, Ml Ml) = - (Ml + M2)2] [p" - (Ml - M2)2' 
For Ml = Ma = M one gets 



(B.4) 



1 



M^ 



f'^ B{M\M\pI = -J{pI + _- In — + 1 



(B.5) 



where J {pi is given by 

j{pI - 



(47r 



1 r /(T + i 



a 



'1- 



4M2 



(B.6) 



The one-loop amphtudes contain an additional logarithmic dependence on the chiral 
renormalization scale u, through the factors (P = it, K, t]): 



Ml 



Up 



In 



Ml 



[B.7) 



C One-loop amplitudes 



The one-loop K ^ rni amplitudes have been computed in refs. |^2|,|4^, in the absence of 
electroweak corrections. The electromagnetic contributions have been recently calculated 



in refs. [^5|-^. The results take the form: 



Al4'^ = i^-\{2Ml-M'^)B{MlMlMl)-^M',B{MlMlMl) 

+ \^{mI- AmI) B{Ml, Ml Ml + j^MkB(Ml, Ml Ml 



4 M2 



+ 



4(Mi - M2)M2 



M2 M^ \ 



Ml ( Ml\ ( 



K ^"^K / 

Ml . Mt 



■Ml 



Ml, 



Ml 



(C.l) 



i (2Mi - Ml) B{Ml Ml Ml) + ^-^B{Ml, Ml, Ml) 



+ - 



IMl 



{Ml - AMI) B{Ml, Ml Ml + ^,B{Ml, Ml Ml 



A Ml 



Ml' 



1/^ Ml\ 1/ Ml\ 3/ 



M]^ 
Ml 



V ( ' 



(C.2) 



for the octet isoscalar amplitude, 



AlAo 



(27) 



-i {2Ml - Ml) B{Ml Ml Ml) + \mIB{MI Ml M\ 
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+ 



+ 



4M2 



{Ml - AMI) B{Ml, Ml Ml) - ]-^B{MI, M^ MI) 



3M2 



Mi 



4(Mi - Ml)Ml 
' Ml 



Ml Mt ' 



10 + 4 



Mi 
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Mi, 



( M2 M^ ^ 



(C,3) 



for the 27-plet isoscalar amplitude and 



\[Ml-2Ml)B{MlMlMl) 



2 

+ - 

+ 



5M2 



K 



Ml 



Ml - ^-Ml) BiMl, Ml Ml + l^^BiMl, M^ Ml 



Mi 



4(Mi - Ml)Ml 



Ml Mt 
5-18^ + 21 



Ml 



Ml, 



,Mr 

'Ml, 



(C.4) 



-(Ml-2Ml)B{MlMlM\ 



K) 



5 Ml 
^ 8 Ml 



Mt 



Ml - -Ml ) B{Ml, Ml Ml + ^B{Ml, M^ Ml 



5 fMl 11 



^ 4 \ M2 



'Ml 



2 \ Ml 



\ 3 fM% 



1 



V ( ' 



(C.5) 



for the 1 = 2 amphtude. 



D Resummation of higher— order corrections 

In this appendix we provide some details on the Omnes procedure for calculating the 
isospin amplitudes. The resummed loop corrections are contained in the factors Cj^\ as 
defined in eq. ([561). At 0{p'^) in the chiral expansion, these quantities should reproduce the 
one-loop xPT results in (^if); this determines the factors C\^\sq), with sq the subtraction 
point, up to higher-order local contributions: 



Cf^(so) = cf^[l-5fi,(Mi,so)J ^ l + ALAr-5niiMl,so). (D.l) 

Here, AlA\^^ is the one-loop xPT result and SQj{Ml, sq) is obtained by taking the chiral 
prediction for the phase shift 6q{z) in Qj{Ml, sq) and expanding Qi{Ml, sq) to 0{pl, 

ni{Ml,so) = l + 6Qj{Ml,so) + 0{p'). (D.2) 
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The explicit expressions for A^.^^^^ are listed in appendix y. The once-subtracted Omnes 
exponential 



n/(Mf,so) = exp 



{Ml - So) 



IT 



dz 



4M2 {z - So) {z- Ml- ie) 



(D.3) 



contains the integral over the experimentally determined phase shift Sq{z) with J = or 2. 
The upper edge of the integral z should correspond to the first inelastic threshold in the 
given isospin channel. The corresponding expansion factor, 



5ni{Ml,so) 



K ~ ■^Oj 



dz 



4M| - So) {z 



Ml 



le) 



(D.4) 



contains the same dispersive integral, but with the phase shift ^0(2;) determined at 0{p^) 
in xPT. The explicit expressions for 5f2/(M|^, sq) with / = and 2 are as follows: 



5no{Ml,so) 



2 Ml -Ml) a{Ml) In 



-{2 So -Ml) a{so) In 



a{Ml)-a{z) 
a{Ml) + a{z) 



aiz 



Sn2iMl,so) 



2 Ml -Ml) a{Ml) In 



a{so) + cr{z 

a{Ml) - a{z) 



2 [Ml -So) In 



(T{Z 



l + cr(z) 



(D.5) 



32 7r2/2 
-{2 Ml -So) a{so) In 



aiMl) + aiz) 



a{so) - a{z) 
c^iso) + cr{z) 



+ [Ml - So) In 



1 - ^(^) 
l + a(z) 



where for convenience we have defined a{s) = yl — AM^/s. 

In the following numerical analysis we have varied the subtraction point between sq 







and So = 3M^, together with the upper edge of the Omnes integral z, to estimate the 
sensitivity of our predictions to these parameters. We have fixed the xPT renormalization 
scale at z/ = Mp. In tables ^ and |^ the dispersive part of the Omnes factors and 6Qj{Ml, sq) 
are reported as functions of Sq, for z = 1 GeV^ and z = 2 GeV^ respectively. The 
corresponding moduli of the corrections C) , derived according to eq. (p^D, are given in 
tables ^ and |. The residual tiny dependence of on the subtraction point sq should 

be cancelled by missing 0{p^) contributions to Cj^\so), since the local ambiguity of the 
Omnes procedure has been only solved to in the chiral expansion. From tables ^ 

and ^ one can also verify that the once-subtracted result is sensitively dependent on z. 

As it was noticed in ref. p[, the sensitivity to the higher energy region of the dispersive 
integral (i.e. the numerical dependence on the upper edge z) is reduced by performing more 
subtractions. However, a better knowledge of the theory is required in this case. Indeed, 
the sensitivity to unknown higher-order corrections in the chiral expansion will increase 
with the number of subtractions, so that the resulting amplitudes can only be trusted at 
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Table 2: The sq dependence of the once-subtracted Omnes parameters for 2; = 1 GeV^. 



So 


% 


6Q0 


3?2 


6Q2 





1.45 


0.32 + 0.46 i 


0.94 


-0.16 - 0.20 i 


Ml 


1.40 


0.29 + 0.46 i 


0.95 


-0.15 - 0.20 i 


TV 


1.33 


0.25 + 0.46 i 


0.96 


-0.13 - 0.20 i 


3Ml 


1.26 


0.21 + 0.46 i 


0.97 


-0.12 - 0.20 i 



Table 3: The Sq dependence of the once-subtracted Omnes parameters for z = 2 GeV^. 



So 


3^0 


5^0 


3^2 







1.58 


0.47 + 0.46 i 


0.92 


-0.24 - 0.20 i 


Ml 


1.51 


0.43 + 0.46 i 


0.93 


-0.22 - 0.20 i 


2Ml 


1.44 


0.39 + 0.46 i 


0.94 


-0.20 - 0.20 i 


ml 


1.35 


0.33 + 0.46 i 


0.95 


-0.17-0.20 i 



Table 4: Resummed loop corrections with one subtraction and z = 1 GeV^. 



So 




1/^(27), 
1^0 \ 


|^(eui)| 
I'-'O 1 


1/^(27), 
l'-'2 1 


l'-'2 1 





1.37 


2.47 


1.38 


1.06 


0.62 


Ml 


1.36 


2.42 


1.37 


1.05 


0.61 


2Ml 


1.35 


2.36 


1.36 


1.05 


0.60 


3Ml 


1.33 


2.28 


1.34 


1.04 


0.59 



the lowest values of the subtraction point (sq ~ 0), where xPT corrections are moderate. 
We have checked these statements using the twice-subtracted Omnes exponential 0]: 

^i{Mk, So) = exp <^ - Sq) — ^ + / 7 r^-, rj^ > , 

[ l + gj{so) IT JiM^ {z - SoY {z - Mf^ - te) } 

(D.6) 

where the functions gi{s) (and their first derivatives g'lis)) are the one-loop contributions 



28 



Table 5: Resummed loop corrections with one subtraction and z = 2 GeV^. 



■So 


I'-'O 1 


I'-'O 1 


I'-'O 1 


|/-(27)| 


\ 





1.26 


2.45 


1.27 


1.10 


0.68 


Ml 


1.27 


2.41 


1.27 


1.10 


0.67 


TT 


1.27 


2.36 


1.28 


1.09 


0.65 


3M2 


1.26 


2.28 


1.27 


1.08 


0.64 



Table 6: The z dependence of the twice-subtracted Omnes parameters for sq = 0. 



z (GeV2) 


3^0 


5^0 


3^2 




1 


1.44 


0.40 + 0.46 i 


0.86 


-0.20 - 0.20 i 


2 


1.46 


0.42 + 0.46 i 


0.85 


-0.21 - 0.20 i 



Table 7: Resummed loop corrections with two subtractions and sq = 0. 



z (GeV2) 


1/^(8), 
I'-'O 1 


1/^(27), 
I'-'O 1 


I'-'O 1 


1/^(27), 
l'-'2 1 


|^(e«))| 
l'-'2 1 


1 


1.25 


2.34 


1.26 


1.00 


0.60 


2 


1.23 


2.34 


1.24 


1.00 


0.61 



(and their derivatives) to the isospin amplitudes due to the elastic vrvr rescattering: 

go{s) = -\{2s-Ml)B{MlMls), 

92{s) = \{s-2Ml)B{MlMls). (D.7) 
The expansion of VLji^Mj^^ sq) at OijP') defines 

5n,{Ml, so) = {Ml - so) g'Aso) + ['^ ... > m) 



TT 



4M2 {z - so)2 {z- Mj.- ie) 



where the phase shift Sq{z) is taken at 0{p'^) in xPT. The numerical results obtained at 
So = 0, with the twice-subtracted Omnes procedure, are reported in tables ^ and |^. 
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The final results for the moduli of the correction factors Cj^\ quoted in eq. (|62D, take 
into account the sensitivity to Sq and z of the once-subtracted solution and the values 
obtained at sq = with two subtractions. 

References 

E. Pallante and A. Pich, Phys. Rev. Lett. 84 (2000) 2568. 
E. Pallante and A. Pich, Nucl. Phys. B592 (2000) 294. 

A. Ceccucci, New measurement of direct CP violation in two pion decays of neu- 
tral kaons by experiment NA48 at CERN, CERN Particle Physics Seminar (February 
29, 2000), |http://www.cern.ch/lNA48/ Welcome.html ; NA48 collaboration (V. Fanti 



et a/.), Phys. Lett. B465 (1999) 335. 

KTeV collaboration (A. Alavi-Harati et al.) Phys. Rev. Lett. 83 (1999) 22. 

A.J. Buras et al., Nucl. Phys. B592 (2001) 55; S. Bosch et al., Nucl. Phys. B565 (2000) 
3; G. Buchalla et al.. Rev. Mod. Phys. 68 (1996) 1125; A.J. Buras et al, Nucl. Phys. 
B408 (1993) 209, B400 (1993) 37, 75, B370 (1992) 69. 



M. Ciuchini et al., Nucl. Phys. B (Proc. Suppl.) 99 (2001) 27; |hep-ph/9910237| ; Nucl. 



Phys. B523 (1998) 501, B415 (1994) 403; Phys. Lett. B301 (1993) 263. 

S. Bertolini et al.. Rev. Mod. Phys. 72 (2000) 65; Phys. Rev. D63 (2001) 056009; Nucl. 
Phys. B449 (1995) 197, B476 (1996) 225, B514 (1998) 63, 93. 

M. Fabbrichesi, Phys. Rev. D62 (2000) 097902. 

T. Hambye et al., |hep-ph/0001088| ; Nucl. Phys. B564 (2000) 391; Eur. Phys. J. CIO 
(1999) 271; Phys. Rev. D58 (1998) 014017; Y.-L. Wu, |hep-ph/001237T| . 

J. Bijnens and J. Prades, JHEP 06 (2000) 035; Nucl. Phys. (Proc. Suppl.) 96 (2001) 
354. 

A. Pich, Nucl. Phys. B (Proc. Suppl.) 93 (2001) 253; E. Pallante, A. Pich and I. 
Scimemi, |hep-ph/00l022g| ; |hep-ph/00l0073| . 



A. A. Belkov et al, |hep-ph/9907335| . 
H.Y. Chen, Chin. J. Phys. 38 (2000) 1044. 
S. Narison, Nucl. Phys. B593 (2001) 3. 

There is a vast literature on this issue, which can be traced back from the recent 
review by A. Masiero and O. Vives, [hep-ph/0104027. 



30 



[16] G. 't Hooft, Nucl. Phys. B72 (1974) 461; B75 (1974) 461. 

[17] E. Witten, Nucl. Phys. B149 (1979) 285; B160 (1979) 57; Ann. Phys. 128 (1980) 363. 

[18] W.A. Bardeen, A.J. Buras and J.-M. Gerard, Nucl. Phys. B293 (1987) 787; Phys. Lett. 
B211 (1988) 343, B192 (1987) 138, B180 (1986) 133; A.J. Buras and J.-M. Gerard, 
Nucl. Phys. B264 (1986) 371. 

[19] K.G. Wilson Phys. Rev. 179 (1969) 1499; W. Zimmermann, Lectures on Elementary 
Particles and Quantum Field Theory, Brandeis Summer Institute (1970) Vol. 1 (MIT 
Press, Cambridge, MA, 1970). 

[20] E.C.G. Stueckelberg and A. Peterman, Helv. Phys. Acta 26 (1953) 499; M. Gell-Mann 
and F.E. Low, Phys. Rev. 95 (1954) 1300; C. Callan Jr., Phys. Rev. D2 (1970) 1541; 
K. Symanzik, Commun. Math. Phys. 18 (1970) 227; 23 (1971) 49. 

[21] A. Pich, Effective Field Theory, in "Probing the Standard Model of Particle Inter- 
actions", Proc. 1997 Les Houches Summer School, eds. R. Gupta et al. (Elsevier, 
Amsterdam, 1999), Vol. II, p. 949. 

[22] M.K. Gaillard and B.W. Lee, Phys. Rev. Lett. 33 (1974) 108; G. Altarelh and L. Ma- 
iani, Phys. Lett. B52 (1974) 351. 

[23] A.I. Vainshtein, V.I. Zakharov and M.A. Shifman, JETP Lett. 22 (1975) 55; Nucl. 
Phys. B120 (1977) 316. 

[24] F.J. Gilman and M.B. Wise, Phys. Rev. D20 (1979) 2392; Phys. Rev. D21 (1980) 
3150. 

[25] A.J. Buras, Weak Hamiltonian, CP Violation and Rare Decays, in "Probing the Stan- 
dard Model of Particle Interactions", Proc. 1997 Les Houches Summer School, eds. 
R. Gupta et al. (Elsevier, Amsterdam, 1999), Vol. I, p. 281. 

[26] A.J. Buras, M. Jamin and M.E. Lautenbacher, Nucl. Phys. B408 (1993) 209; 
Phys. Lett. B389 (1996) 749. 

[27] M. Ciuchini et al, Phys. Lett. B301 (1993) 263; Z. Phys. C68 (1995) 239. 

[28] S. Weinberg, Physica 96A (1979) 327. 

[29] J. Gasser and H. Leutwyler, Ann. Phys., NY 158 (1984) 142; Nucl. Phys. B250 (1985) 
456; 517; 539. 

[30] G. Ecker, Prog. Part. Nucl. Phys. 35 (1995) 1. 

[31] A. Pich, Rep. Prog. Phys. 58 (1995) 563. 

[32] U.-G. Meissner, Rep. Prog. Phys. 56 (1993) 903. 



31 



[33] J. Gasser, U-G. Meissner, Phys. Lett. B258 (1991) 219. 

[34] J.F. Donoghue et al., Phys. Lett. B179 (1986) 361; H.Y. Cheng, Phys. Lett. B201 
(1988) 155; M. Lusignoh, Nucl. Phys. B325 (1989) 33. 

[35] A.J. Buras and J.-M. Gerard, Phys. Lett. B192 (1987) 156. 

[36] G. Ecker, G. Miiller, H. Neufeld and A. Pich, Phys. Lett. B477 (2000) 88. 

[37] J. A. Cronin, Phys. Rev. 161 (1967) 1483. 

[38] J. Bijnens and M.B. Wise, Phys. Lett. B137 (1984) 245. 

[39] B. Grinstein, S.-J. Rey and M.B. Wise, Phys. Rev. D33 (1986) 1495. 

[40] C. Bernard et al., Phys. Rev. D32 (1985) 2343. 

[41] R.J. Crewther, Nucl. Phys. B264 (1986) 277. 

[42] J. Kambor, J. Missimer and D. Wyler, Nucl. Phys. B346 (1990) 17; Phys. Lett. B261 
(1991) 496; J. Kambor et al, Phys. Rev. Lett. 68 (1992) 1818. 

[43] J. Bijnens, E. Pallante and J. Prades, Nucl. Phys. B521 (1998) 305; E. Pallante, JHEP 
01 (1999) 012. 

[44] A. Pich, B. Guberina and E. de Rafael, Nucl. Phys. B277 (1986) 197. 

[45] G. Ecker, G. Isidori, G. Miiller, H. Neufeld and A. Pich, Nucl. Phys. B591 (2000) 419; 
work in progress. 

[46] V. Cirigliano, J.F. Donoghue and E. Golowich, Phys. Lett. B450 (1999) 241; Phys. 
Rev. D61 (2000) 093001, 093002. 

[47] V. Cirighano and E. Golowich, Phys. Lett. B475 (2000) 351. 

[48] G. Ecker, J. Kambor and D. Wyler, Nucl. Phys. B394 (1993) 101. 

[49] E. de Rafael, Nucl. Phys. B (Proc. Suppl.) 7A (1989) 1. 

[50] A. Pich, Nucl. Phys. B (Proc. Suppl.) 7A (1989) 194. 

[51] A. Pich and E. de Rafael, Nucl. Phys. B358 (1991) 311; Phys. Lett. 374 (1996) 186. 
[52] M. Jamin and A. Pich, Nucl. Phys. B425 (1994) 15. 

[53] G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B321 (1989) 311; G. Ecker, 
J. Gasser, H. Leutwyler, A. Pich and E. de Rafael, Phys. Lett. B223 (1989) 425. 

[54] S. Peris, M. Perrotct and E. de Rafael, JHEP 05 (1998) Oil. 



32 



M. Jamin, J. A. Oiler and A. Picli, Nucl. Phys. B587 (2000) 331; work in progress. 

G. Amoros, J. Bijnens and P. Talavera, Phys. Lett. B480 (2000) 71; Nucl. Phys. B585 
(2000) 293. 

J. Bijnens, G. Colangelo, G. Ecker, JHEP 02 (1999) 020. 

M. Davier, S. Chen, A. Hocker, J. Prades and A. Pich, Strange quark mass from t 
decays, Proc. 6th Int. Workshop on Tau Lepton Physics (Victoria, Canada, 2000), 
Nucl. Phys. B (Proc. Suppl.) 98 (2001). 

A. Pich and J. Prades, JHEP 10 (1999) 004, JHEP 06 (1998) 013. 
ALEPH Collab., Eur. Phys. J. Cll (1999) 599, CIO (1999) 1. 

J. Kambor and K. Maltman, Phys. Rev. D62 (2000) 093023; K.G. Chetyrkin et al., 
Nucl. Phys. B533 (1998) 473. 

M. Jamin, Nucl. Phys. B (Proc. Suppl.) 64 (1998) 250; K.G. Chetyrkin et al., 
Phys. Lett. B404 (1997) 337; P. Colangelo et al, Phys. Lett. B408 (1997) 340; 
C.A. Dominguez et al., Nucl. Phys. B (Proc. Suppl.) 74 (1999) 313; K. Maltman, 
Phys. Lett. B462 (1999) 195; S. Narison, Phys. Lett. B466 (1999) 345. 

V. Lubicz, Nucl. Phys. B (Proc. Suppl.) 94 (2001) 116; R. Gupta and K. Maltman, 
|hep-ph/0 10113^ . 

J. Prades, Nucl. Phys. B (Proc. Suppl.) 64 (1998) 253; J. Bijnens, J. Prades and E. de 
Rafael, Phys. Lett. B348 (1995) 226. 

A. Pich, Tau Physics: Theoretical Perspective, Proc. 6th Int. Workshop on Tau Lep- 
ton Physics (Victoria, Canada, 2000), Nucl. Phys. B (Proc. Suppl.) 98 (2001) 385; 
Tau Physics, Proc. XIX International Symposium on Lepton and Photon Interactions 
at High Energies (Stanford, 1999), eds. J. Jaros and M. Peskin (World Scientific, 
Singapore, 2000) 157 ||hep-ph/991229l . 



Particle Data Group, Review of Particle Physics, Eur. Phys. J. C15 (2000) 1; and 
2001 update ||http://pdg.lbl.gov/ 



G. Rodrigo, A. Pich and A. Santamaria, Phys. Lett. B424 (1998) 367. 

K.G. Chetyrkin, B.A. Kniehl and M. Steinhauser, Phys. Rev. Lett. 79 (1997) 2184; 
W. Bernreuther and W. Wetzel, Nucl. Phys. B197 (1982) 228; W. Wetzel, Nucl. Phys. 
B196 (1982) 259; W. Bernreuther, Ann. Phys. 151 (1983) 127. 

K. Maltman and C.E. Wolfe, Phys. Lett. B482 (2000) 77; Phys. Rev. D63 (2001) 
014008. 

G. Ecker, A. Pich and E. de Rafael, Phys. Lett. B237 (1990) 481. 



33 



71] S. Gardner and G. Valencia, Phys. Lett. B466 (1999) 355. 

72] F. Guerrero and A. Pich, Phys. Lett. B412 (1997) 382. 

73] R. Omnes, Nuovo Cimento 8 (1958); 

N.L Muskhelishvili, Singular Integral Equations, Noordhoof, Groningen, 1953. 

74] T.N. Truong, Phys. Lett. B207 (1988) 495. 

75] M. Biichler et al, |hep-ph/010228^ , |hep-ph/010228£ . 

76] M. Ciuchini et al, NucL Phys. B573 (2000) 201, |hep-ph/0012308| ; S. Mele, Phys. Rev. 
D59 (1999) 113011; F. Parodi, P. Roudeau and A. Stocchi, Nuovo Cim. A112 (1999) 
833; A. Hocker, H. Lacker, S. Laplace, F. Le Diberder, [hep-ph/0104062 . 

77] J. Prades and A. Pich, Phys. Lett. B346 (1995) 342; J.L. Rosner, |hep-ph/001lT8^ . 

78] A. Buras, M. Jamin and P.H. Weisz, Nucl. Phys. B347 (1990) 491; S. Herrlich and 
U. Nierste, Nucl. Phys. B419 (1994) 292, B476 (1996) 27; Phys. Rev. D52 (1995) 6505. 

79] L. Wolfenstein, Phys. Rev. Lett. 51 (1983) 1945. 

80] E. de Rafael, Chiral Lagrangians and Kaon CP-Violation, in "CP Violation and the 
limits of the Standard Model" (TASr94), ed. J.F. Donoghue (World Scientific, Singa- 
pore, 1995) ||hep-ph/950225l . 

81] S. Peris and E. de Rafael, Phys. Lett. B490 (2000) 213. 

82] J. Bijnens, H. Sonoda and M.B. Wise, Phys. Rev. Lett. 53 (1984) 2367. 

83] M. Knecht, S. Peris and E. de Rafael, |hep-ph/01020T7| ; Nucl. Phys. B (Proc. Suppl.) 
B86 (2000) 279; Phys. Lett. B457 (1999) 227. B443 (1998) 255. 



84] L. Lellouch and M. Liischer, |hep-lat/0003023| ; M.F.L. Golterman and E. Pallante, 
JHEP 08 (2000) 023; C.-J.D. Lin, G. Martinelh, C.T. Sachrajda and M. Testa, f^- 
lat/0104006| . 



34 



